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ABSTRACT 

We present Spitzer MIR spectra of 25 FR~I radio galaxies and investigate the nature of their 
MIR continuum emission. MIR spectra of star-forming galaxies and quiescent elliptical galaxies 
are used to identify host galaxy contributions while radio/optical core data are used to isolate the 
nuclear non-thermal emission. Out of the 15 sources with detected optical compact cores, four 
sources are dominated by emission related to the host galaxy. Another four sources show signs 
of warm, nuclear dust emission: 3C15, 3C84, 3C270, and NGC6251. It is likley that these warm 
dust sources result from hidden AGN of optical spectral type 1. The MIR spectra of seven sources 
are dominated by synchrotron emission, with no significant component of nuclear dust emission. 
In parabolic SED fits of the non-thermal cores FR-Is tend to have lower peak frequencies and 
stronger curvature than blazars. This is roughly consistent with the common picture in which 
the core emission in FR-Is is less strongly beamed than in blazars. 

Subject headings: galaxies: active — galaxies: jets — galaxies: nuclei 



1. Introduction 

Large double radio sources come in two radio- 
morphological types. The more powerful type 
is the "classical double" or FR-II source, whose 
prototype is Cygnus A. This group has strong 
terminal "hotspots" , and generally one-sided 
arcsecond (kpc) scale. The other type, FR- 
Is like Centaurus A and M87, is lower in ra- 
dio luminosity, double-j etted, and edge-darkened 
( Fanaroff &: RilevI Il974l ). These sources mostly 



have symmetric kpc jets. 

It has been known since the 1960s that some of 
these giant (> 100 kpc) radio sources show pow- 
erful optical/UV continuum emission (called "Big 
Blue Bump" , or BBB) from a central point source. 
This optical/UV emission is generally identified 
as thermal radiation from dense matter accreting 
onto supermassive black holes. The gravitational 
potential energy of this infalling matter is thought 
to power much or all of the observed activity. Ob- 
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jects showing this strong optical/UV continuum 
and accompanying broad emission Hnes are classi- 
fied as quasars and broad Une radio galaxies (here- 
after, "quasars"). On the other hand, double ra- 
dio sources without these components are called 
(Narrow Line) Radio Galaxies (NLRGs). 

Many FR~II NLRGs have hi dden quasars de- 
tected by spectropolarimetry (e.g. Antonuccil 19841 : 



Ogle et al.Hl997t ICohen et al.lll999l) and recently 



the mid infrared (e.g. 


Meisenheimer et al. 


2001 


Haas et aL,,2004;iS 


bi et al.i,2005: Ode et al. 


2006 


Clearv et all 2007r 


Tadhunter et al. 


|2007|). 


Thus 



many powerful radio galaxies contain quasars 
which are hidden from direct view by dusty equa- 
torial tori - i.e . as expected in the "Unified Model" 
(iBarthell [19891 : lAntonuccilll993t luTrv fc Padovani 



19951) 



The total light optical spectra of many FR-II 
NLRGs already indicate the presence of a "classi- 
cal" AGN (those thought to accrete radiatively ef- 
ficiently) by their high ionization emission line ra- 
tios (high ionization galaxies - HIG(1]). Some FR- 
II RGs and most FR-Is, however, show only weak 
line emission of low ionization and are therefore 
classified as low ionization galaxies (LIGs). Here 
a luminous AGN (if present) has to be discovered 
at other wavelengths. Alternatively, these radio 
galaxies may simply lack a hidden quasar. A lack 
of a BBB means, according to current wisdom, a 
lack of a copious and radiatively efficient accretion 
fiow. Historically such a situation was thought 
to require energy extraction from black hole rota- 
tion, and var ious mechanisms have been suggeste d 
for this (e.g. iBegelman. Blandford. fc Reed 1984[ ). 



Alternatively in principle accretion power could 
make the jet, yet not radiate significantly. In such 
a case the BBB is present, but of low power com- 
pared to a (lobe-matched) quasar BBB. 

In an optical imaging study o f 3CR FR-I 



minosity is taken as evidence for their common 
origin as synchrotron emission from the base of a 
jet. Assuming that the base of the jet is found on 
scales smaller than any obscuring torus and con- 
sidering the high detection rate of CCCs (85 %) 
they argue that we see directly into the nucleus 
and thus no standard, geometrically thick torus 
can be present in most low-luminosity (FR-I) 
radio galaxies. Since no broad lines or strong 
BBB emission is observed there could also be no 
hidden, radia tively efficient AGN in those cases. 
However, e.g. Cao fc Rawlingsl ( 2004h note that if 
the emission from the jet occurs on scales larger 
than that of the torus the absence (or presence) 
of a torus cannot be directly inferred from the 
observed nuclear jet components. In addition, 
high-resolution CO maps have revealed massive 
nuclear tori like those invoked in the standard 
unified model in the ce nter of at least two FR- 
I radio gal axies: 3C31 (lOkuda et~aL 2005 ) and 
NGC 3801 (|Das et al.ll2005l ). Thesetori have high 
molecular masses (>10^ M©), with high molecular 
column densities and probably geometrically thick 
shapes. 

From an UV imaging study IZirbel fc Baum 



(|2003[ ) argue that the properties of FR-Is are con- 
sistent with the classic unification models and the 
existence of an obscuring torus: While BL Lac ob- 
jects are well aligned sources, FR-I sources with 
nuclear UV components are at a critical angle for 
which a torus hides the nucleus but not the jet 
(or the latter greatly outshines the former). FR- 
Is without nuclear UV components arc thought to 
have a larger viewing angle where a torus hides 
the central regions completely. 

Recently, some of the compact optical cores dis- 
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covered bv lChiaberge et aL ( 19991 ) have been stud 
led with imaging polarimetry revealing fairly large 
polarization (< 11%) which further indicates their 
non-thermal origin, though such high polariza- 



(Il999l) report unresolved central compact cores 



source s using Hubble Space Telescope, IChiaberge et all tions can also result from scattering ([Capetti et al 



20071) 



(CCCs) in the centers of most observed galaxies. 
The correlation of optical (at ~ 7000 A) compact 
cores and radio (at 6 cm) cores in flux and lu- 



^ Although such objects are usually referred to as high exci- 
tation galaxies (HEGs) the term HIGs seems more appro- 
priate because the distinction is primarily one of ionization 
and not excitation, as expected for photoionized gases. 



In the X-rays, FR-I sources are dominated 
by weakly absorbed non-thermal jet compo- 
nents in most case s (e.g. iBalmaverde et al. I l2006hl: 
lEvans et"al] l2006f) . They lack the powerful ab- 
sorbed, accretion-related component of a strong 
AGN (associated with a luminous accretion disk 
and circumnuclear obscuring structure). In many 
FR-IIs the latter is found in addition to the 
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non-thermal jet component (jEvans et al. l2006h . 
Thus, in FR-Is any non-jet accretion-related X- 
ray component carries much lower intrinsic AGN 
luminosities (10^^-10^^ ergs~^) or else their detec- 
tion is hampered by column densities of well above 
10^4 cm-2. Also IWu et all (|2007l ) argue that the 
(soft) X-rays in their FR-I sample can be ac- 
counted for by jet emission only. We must note 
however that FR-Is with broad emission lines in 



the o ptical do e xist (e.g. 3C120. iTadhunter et al 
1993 : see also AntonuccH l2002al for an anec- 



dotal listing of such objects). Such broad- 
line FR-Is can extend up t o quasar luminosities 
( Gower k, Hutchind Il984t iBlundell fc Rawlingd 



2001 



Hevwood et"aLH2007l ). 



We here present a study of the mid-infrared 
properties of a sample of FR-I radio galaxies. 
The motivation is to possibly detect a hidden 
AGN by the reprocessed warm thermal dust emis- 
sion as found in many FR-II radio galaxies (e.g. 
Ogle et ahl l2006l ). And even in the case of ex- 



ceptionally high extinction which attenuates also 
the MIR luminosity, spectral features that give 
clues on the thermal or non-thermal origin of 
the emission might still be detectable. On the 
other hand, the nuclear part of the MIR spectra 
bridges the gap in wavelength between the opti- 
cal compact cores and the radio cores observed by 
IChiaber ge et al. ( 1999[ ). If there is in fact no hid- 
den AGN and no obscuring structure to be found 
in FR-Is, the MIR spectra should agree with syn- 
chrotron core estimates. 

2. Sample selection and observations 

Our sample of FR-I radio galax ies was selected 
from the 3CR ( Spinrad et al.lll985 ) catalog (about 
half of the sources are also included in the SCRrEI). 
We selected a tractable subset of lobe dominated 
sources with a flux limit of 15 Jy at 178 MHz. To 
increase the sample size we added some very sim- 
ilar sources which were close to, but technically 
below, the flux cut off. We also observed two ad- 
ditional sources because they were of special inter- 
est: 3C18d!l as well as IC 4296. See Tab.[I]for some 
basic data of the selected sources. We note that 



^ http: / / 3crr . extragalactic. info / 

33C189 is only part of the 3C catalog llEdee et al.lll959f) 
but not of the S CR or the 3CRR catalogues. See e.g. 
iLaing et al] ||1983| ') for details on these samples. 



our sample is not complete but can be considered 
representative for lobe dominated FR-I sources. 
No sources were lost by imposing scientifically rel- 
evant selections. 

The morphological classification of the sources 



accor ding to the Fanaroff-Riley scheme (Fanar off fc Rilev 
19741 ) was checked using radio maps from the lit- 
erature (see Tab. 1 for references). 

The objects w ere observed in low -resolution 



mode with IRS (jHouck et al.l 120041) on boar d 
Spitzer Space Telescope (| Werner et al.l 120041 ) . 



mostly under pid-20525. Each source under this 
pid was observed in staring mode for a total of 960 
sec, 1440 sec, 720 sec, and 480 sec in SL2, SLl, 
LL2, and LLl, respectively. After averaging all 
the cycles at one nod position we subtracted the 
off-order frames to remove the background. The 
resulting images were cleaned of residual rogue 
pixels and cosmic rays using Irsclean. The 
Spice software was used to calibrate and extract 
one dimensional spectra. The extracted spectra 
from the two nod positions were then averaged 
and data from the different modules were com- 
bined into a single spectrum. Due to the larger 
slit width in LL (10.7" in LL versus 3.7" in SL), 
multiplicative factors (< 1.5) were applied to the 
SL modules where necessary in order to match the 
flux levels. 

Spectra for additional FR-I sources not ob- 
served under pid-20525 (see Tab.[Tl lower part) 
were also retrieved from the archive and were sub- 
ject to the same procedures. 

3. The Spitzer spectra 

In FiguresH] and [5] we present the rest frame 
MIR spectra for our FR-I core sample and the 
supplemental objects. Because for 3C15 and 3C29 
the SL modules of IRS suffer from saturated peak- 
up areas we present their spectra separately in 
Fig.[3l Statements on the spectral behavior of the 
individual sources refer to the figures where their 
observed spectra are shown (Figs.[T] and [2]) and 
are thus considered in terms of F^ . All sources in 
our sample have been detected except for 3C403.1 
which was not detected in either SL or LL. 

3.1. The main sample 

3C31 (NGC383) This source shows very 
strong PAH emission and the MIR spectrum is 
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Table 1: The Sample. 



Object 


alt. 


z 


Fl78MHz 


total 


F5 GHz 


core 


CCC 




map 


spec 




name 




Jy 


reP 


mJy 


reP 




filter 


reP 


reP 


3C15 




0.073000 


17.2 


42 


32 


25 


38 


F160W 


22 


48 


3C29 




0.045031 


16.5 


42 


93 


24 


10 


F702W 


24 


48 


3C31 


NGC383 


0.017005 


18.3 


43 


92 


12 


29 


F702W 


19 


48 


3C66B 




0.021258 


26.8 


43 


182 


12 


97 


F702W 


14 


48 


3C76.1 




0.032489 


13.3 


43 


10 


37 






21 


48 


3C83.1 


NGC 1265 


0.025137 


29.0 


43 


25 


27 


3 


F702W 


27 


48 


3C129 




0.020800 


51.1 


44 


34 


38 






36 


48 


3C189 


NGC 2484 


0.040828 


7.4 


45 


195 


13 


76 


F814W 


7 


1 


3C264 


NGG 3862 


0.021718 


28.3 


43 


200 


12 


248 


h 791 W 


41 


28 


3C270 


NGC 4261 


0.007465 


56.5 


42 


308 


24 


11 


F791W 


9 


33 


3C272.1 


M 84 


0.003536 


21.1 


43 


180 


12 


139 


F814W 


18 


48 


IC 4296 




0.012465 


16.8 


46 


297 


24 


524 


F160W 


24 


40 


3C293 




0.045034 


13.8 


43 


100 


12 






5 


34 


3C317 




0.034457 


53.4 


42 


391 


24 


23 


F814W 


6 


48 


NGC6251 




0.024710 


11.6 


43 


350 


16 


191 


F812W 


30 


23 


3C386 




0.016885 


26.1 


43 


14 


32 


2301 


F702W 


21 


48 


3C403.1 




0.055400 


14.7 


42 














3C424 




0.126988 


15.9 


42 


18 


8 






3 


48 


3C465 


NGC 7720 


0.030221 


41.2 


43 


270 


12 


36 


F702W 


10 


11 


3C84 


NGC 1275 


0.017559 


66.8 


43 


3100 


26 


3489 


F702W 


29 


48 


3C120 




0.033010 


7.4 


45 


1970 


39 






31 


33 


3C218 


HydA 


0.054878 


225.7 


46 


217 


24 






35 


40 


3C274 


M87 


0.004360 


1144.5 


43 


4000 


12 


894 


F814W 


15 


48 


BL Lac 




0.068600 


-2.5^= 


47 


-4.0^^ 








2 


20 


E1821+643 




0.297000 


-0.5 


47 


8.0 


4 






4 
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"References for total radio fluxes, radio core fluxes, radio maps, and optical spectra, respectively. 

''Data take n from IChiaberge et al.l lll999l ') except for 3C189 jCapetti et al.l |2002| ). IC4296 llBalmaverde et al.l [2006al ). and 
NGC 6251 l lChiaberge et al.ll2003t) Data on 3C15 provided by R. Baldi (private communication). 

■^The flux of only the extended emission is 29 mJy at 20 cm l lAntonuccilll98d) which correspond to ~ 151 mJy at 178 MHz (using 
a = —0.8; S oc u"). For the variable 5 GHz core flux we here give the average of all measurements available from NED. 



References. — ( 1) SDSS; (2) lAntonuc'cilll986l: fS) iBlack et al.lll992l: (41 iBlundeU fc Rawlingsl [ 200ll: (51 iBridle et al.l[l98ll: 
(6) lBurn slll99(]|: (7) ICapetti et al .lll993l: CSl IChiaberge et al.ni999l: (9')'Condon Sz Brodericklll98'8r(10)ICond on et al.lll99lt (11 



Dc Robcr tis fc YeeHl99(]l: (1 2) [Giovannini ct al. 1988; (13) Giovannini et al. 2005; (14) Hardcastlc ct al. 1996i: (151 iHines et al 
"989 l; (16) 1 Jones et al.lll986l; (17llKol latschnv ct al. 2006; (18) Laing & Bridle 1987; (19) Laing ct al. 2008; (201 iLawre nce et al 



199d: (211 iLeahv fc Perleyl Il99ll: ( 22) Leahy ct al. 1997; (23) Miley & Osterbrock 1979; (24) Morganti ct al. 1991 (25) 
Morganti et al.lll999l ; (261 iTavlor et al. 2006; (27) O'Dca & Owen 1986; ( 281 lOwen et al.lll99Q; (291 iPedlar et al.lll990l: (30 



Morganti et al.lll999l ; (26)_lTavlor et al. 2006; (27) O'Dca & Owen 1986; ( 281 lOwen et al.lll99Q; (29) iPedlar et al.lll990l: (301 
Perlev et al.lll984|: (31 1 [Kmpson ct al. 1996; (32) Strom ct al. 1978; (3 31 iTadhunter et aljll993l ; (34 1 iTadhunter et al.l [ 2005t 
(35) 'Taylor et al."l990'; (36) Taylor ct al. 2001; (37) Valloc 1982; (38) van Brcugcl fc Jagerslll982l; (391 IWalker et al.lll987l; (401 



Wills ct al. 2004; (41) sec littp://www.jb. man. ac.uk/atlas/ (42) KcUcrmann ct al.l | |1969|1 . on lBaars et al.l (1 19771 ) scale using 
Laing fc Pe acock ( 19 8C1|1: (431 taken from http: / /3crr.cxtragalactic.info/ , (44) estimated using GMRT measurements (|Lal fc Rad 



200411 ; (451 Kuchr ct anil98ll ; (46) estimat ed using the 160 MHz flux llKuchr ct al... 1981.1 and a = -0.8; (47) estimated from 



low frequency measurements in NED; (48) iButtiglione et al.ll2009l 



likely to be dominated by emission related to star 
formation. 

A massive molecul ar gas disk lia,s been detected 
in CO in this object (|Okuda et al.ll2005h . 



3C66B This source shows increasing flux to- 
wards the shortest wavelengths observed (starting 
at around 8 fim) indicating a dominant contribu- 
tion from the stellar component of the host galaxy 
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Fig. 1. — Observed MIR spectra of the FR-I sources. Flux in Jy is plotted over the rest frame wavelength 
in fim. 



at these wavelengths. A weak silicate feature in 
emission can be seen around 9.7 /xm as well as 
significant red continuum at longer wavelengths. 
Atomic line emission from low ionization species 
is detected but no PAH features. 

3C76.1 While this source is clearly detected 



(as seen from the 2-D frames) it has low flux and 
shows a noisy spectrum. No spectral features can 
be securely identified. The spectrum shows slowly 
increasing flux towards longer wavelengths, start- 
ing at around 15 /^m. 

3C83.1 (NGC1265) Emission from the 
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Fig. 1. — continued 



host galaxy stellar population dominates this 
source up to ~ 15 /xm where the spectrum starts 
to flatten in flux. PAH emission at 11.3 /zm and 
atomic emission from [Ne II] at 12.8 /im arc clearly 
detected. The ratio between the PAH features at 
7.7 /um and 11.3 ^m appears very small compared 
to e.g. 3C31 (see §4.1). 



3C129 Here the stellar emission again dom- 
inates below ~ 10 fim and for A > 15 fiia strong 
MIR continuum emission is observed. We also de- 
tect weak silicate emission and PAHs at 11.3 /im. 

3C189 (NGC2484) This source appears 
very similar to 3C66B in the MIR spectrum with 
a distinct stellar component at short wavelengths, 
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Fig. 1. — continued 



weak silicate emission around 9.7/im, red contin- 
uum emission for A > 10 /im, and no PAHs. 

3C264 (NGC3862) Strong PAH features 
are detected in this source which argues for a con- 
siderable, if not dominant, contribution from emis- 
sion related to star formation. In fact, the steep 
continuum at A > 20 /im might well be associ- 
ated with star formation. At the shortest wave- 
lengths some minor contributions from the host 
galaxy stars can be identified. 

3C270 (NGC 4261) A weak silicate feature 
and weak PAH emission is detected. While the 
blue spectral slope towards shorter wavelengths in- 
dicates contributions by stellar emission, a strong 
red continuum is observed for A > 10 /im. In com- 
bination with the relatively weak PAH emission 
this continuum emission is not likely to originate 
from star-forming activity only. 

3C272.1 (M84 - NGC 4374) A very un- 
usual spectrum which has a blue slope up to 



^ 20 /im and shows a very steep red slope at 
A > 30 /tm. Strong PAH emission is observed, but 
as in the case of 3C83.1 the 7.7 /im to 11.3 /im ratio 
appears very small. Several atomic emission lines 
are detected, among them [O IV] , a relatively high 
ionization line which can, ho weve r, also be found 
in star-forming galaxies (e.g. ISmith et al.ll2007t l. 

IC 4296 Somewhat similar in appearance to 
3C272.1 the continuum shows a red slope for A > 
20 /im which steepens for A > 30 /im. Weak PAH 
emission and strong atomic line emission is de- 
tected. The [O IV] emission line has much larger 
equivalent width than in 3C272.1 although the 
relative contributions from star formation to the 
spectrum, as traced by the PAHs, is apparently 
lower. 

3C293 The MIR spectrum of this object is 
clearly dominated by star formation. The strong 
PAH features and the shape of the red contin- 
uum are very similar to what is seen in local star- 
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Fig. 2. — Observed MIR spectra of the supplemental objects. Again we plot the flux in Jy over the rest 
frame wavelength in /im. 



forming galaxies (e.g. lSmith et al ] l2007h . Notably, 
we also detect several rotational transitions from 
molecular hydrogen in emission (Ogle et al., in 
preparation) . 

3C317 This source shows a blue slope at the 
lowest wavelengths suggestive of stellar emission. 



Starting at around 10 /im the spectrum turns to 
a red, but not very steep slope. We do not de- 
tect significant PAH features and [Ne II] A11.3/im 
is the strongest atomic emission line. We do, how- 
ever, detect faint emission from molecular hydro- 
gen from S(l) up to at least S(6) (Ogle et al., in 
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Fig. 3.— The MIR spectra of 3C15 (left) and 3C29 (right). 



preparation). 

NGC 6251 Strong siHcate emission at 9.7 /xm 
and 18 /irn . is d etected as previously noted in 



Ogle et al.l (|2007l ). Strong atomic and weak PAH 



emission can also be seen. 

3C386 This source shows a very starlight- 
dominated, blue continuum throughout the whole 
covered wavelength range. Some PAH features 
(e.g. 7.7/im and 11.3 /im) and atomic lines (e.g. 
[Ne II] A12.8 fim and [Ne III] A15.6 /im) are present 
as well as weak m olecular hydrogen emission. 
ISimpson et al. ( 1996[ ) tentati vely detected a broad 



Ha line in the optical and Madrid et al. ( 20061 ) 



find that a foreground star falls on top of the 
nucleus for this galaxy (previously also noted by 
Lvndslll97lh . 

3C424 The continuum emission is rather 
weak in this source but we detect it in all or- 
ders. The most notable feature in this spectrum 
is the exceptionally bright S(l) emission line from 
molecular hydrogen. Several other H2 emission 
lines can be identified as well. 

3C465 (NGC 7720) While the presence of 
stellar emission can be identified at short wave- 
lengths, the spectrum shows a prominent red slope 
for A > 20 /Ltm. Some atomic emission lines are de- 
tected on top of the rather noisy continuum and 
no significant PAH features (besides a possible de- 
tection of a weak 11.3 /xm feature) are found. In- 
spection of the 2-D spectral frames confirms that 
the source is detected in all orders, but especially 
in LL2 (~ 15 - 20 ^ni) it is of very low S/N. 



3.2. The supplemental sample 

In Fig. [2] we show the spectra of the FR-I 
sources that do not belong to our initial sample. 
We use them here to demonstrate the diversity of 
MIR spectra seen for FR-I sources. 

3C84 (Per A - NGC 1275) This source is 



highly core dominated in the radio (jPedlar et al 



990) 



,19901) . In the MIR, however, a strong thermal 
spectrum, including weak silicate emission and a 
steepening around 12 /im is observed. It appears 
very similar in shape to the spectrum of the clas- 
sical hidden AGN source 3C405 (=CygA; Fig|4]) 
but the continuum appears much stronger in 3C84 
relative to the emission lines. 

3C120 This source is the only 3C FR-I 
source in our sample where an AGN is already 
clearly visible at optical wavelengths as traced 
by broad emission lines and a blue optical/UV 
spectrum (Tab. 1). The Spitzer spectrum re- 
flects this by showing silicate emission at 9.7 ^m 
and 18 /im, a rather flat continuum compared to 
sources like 3C84 and 3C405 and strong [Ne V] 
and [O IV] emission simila r to the MIR spectra 



ipectra 

of typical type-1 AGN fe.g. lWeedman et allliool 
Buchanan et~alll2006l ). 



3C218 Hydra A is the only source in our 
sample which is morphologically clearly a FR-I 
but with an 178 MHz hmrinosity that places it 
among the FR-IIs. Optically a low ionization 
galaxy (Tab. 1) the MIR spectrum shows PAH 
emission and [O IV] is also detected. Overall 
the spectrum looks quite similar to sources dom- 
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inated by star formation in the MIR (e.g. 3C31, 
3C293). However, here the continuum starts rising 
already at around 10 to 12 ^m while for the other 
star-forming sources such a steepenin g in the con- 
tinuum is observed around ~ 20 /im (jSmith et al 



3C405 - Cyg A 



[2007) . Thus it might have a warmer dust compo- 
nent or stronger non-thermal core contributions. 

3C274 (M 87 - Vir A - NGC 4486) The 
MIR spectrum of M87 shows a rather flat slope 
overlayed with some line emission from generally 
low ionization species (except for some weak [O IV] 
emission). At short wavelengths some contribu- 
tions from the host galaxy can be identified and 
there is even some shallow silicate emission de- 
tected. For this source high-resolution ground 
based MIR flux measure ments at ^0.5'^ reso- 
lution are also available (Perlman et al. 2001 



Whvsong fc Antonucci 200 



BL Lac As is common to sources in the class 
of objects named after this prototype, BLLac is 
dominated by Doppler boosted synchrotron emis- 
sion at least from the radio through the optical. 
The MIR spectrum is also dominated by this syn- 
chrotron emission and shows a virtually feature- 
less power-law continuum (a ^—0.7). The dif- 
fuse ra dio emissioii sugge sts a relatively weak FR-I 
source ( Antonuc"cilll986[ ). We show this object to 
illustrate a purely non-thermal spectrum. 

E1821+643 This source was first classified 
as a radio-quiet QSO, but di scovered in a sur- 
vey at low radio frequencies ( Lacv et all Il992). 
Deep radio imaging, however, revealed a ~ 300 kpc 
FR-I radio structure assoc iated with this source 
(jBlundell fc RawlingsllioOlh . Being a broad-line 
object the MIR spectrum looks quite similar to 
those of most other quasars or broad-line objects 
(e.g. 3C120). We clearly detect silicate emission 
at 9.7 /zm and a weak emission feature at 18 /im. 
We also detect [Ne V] emission and weak PAH 
emission (e.g. at 11.3 /im). The continuum long- 
ward of ~ 10 nm is redder than for 3C120 but not 
as red as for hidden AGN sources like Cyg A. 

In Fig. [3] we present the MIR spectra of 3C15 
and 3C29. The spectra of these objects suffer 
from "spill-over" due to saturated peak-up areas. 
Therefore the flux levels and the slopes in the SL 
modules are likely to be compromised. It is no- 
table, however, that both sources clearly show a 
red continuum at wavelengths greater than 15 /im 
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Fig. 4. — Observed MIR spectrum of 3C405, a 
prototypical hidden AGN source. 



but no significant PAH emission. Since the red 
continuum is detected in the LL orders it is not 
an artifact of the spill-over. 

4. MIR emission in FR I radio galaixies 

In this paper we study the nuclear MIR emis- 
sion of FR-I radio galaxies and try to determine its 
properties and possible origins. Components con- 
sidered are i) synchrotron emission from a non- 
thermal core, 11) warm dust emission resulting 
from AGN heating. III) dust heated by star for- 
mation as traced by PAH features, and iv) stellar 
emission (photospheric emission plus dust emis- 
sion associated with the circum-stellar envelopes 
of AGB stars) . We often summarize contributions 
III and IV as processes from the host galaxy. We 
use the term "nuclear" to address central emission 
related directly to an 

In the framework of models which favor the ab- 
sence of nuclear dust with high covering factors 
in FR-I sources, the MIR spectra should consist 
of a combination of non-thermal emission due to 
a nuclear synchrotron source and more extended 
emission related to processes in the host galaxy. 
The warm/hot dust component related to an AGN 
should be missing or of very low luminosity. In 



*We note that the definite test for the nuclear origin of the 
detected MIR emission would be high spatial resolution 
MIR imaging in order to compare the fluxes of any sub- 
arcs ec nuclear point source s with the Spitzer fluxes (see 
e.g. IWhvsone &: Antonuccil ||2004| ') for a discussion of the 
contrasting results on M 87 and Cen A) . 
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Fig. 5. — 3C270 as an example for the subtrac- 
tion of the star-forming template (thin solid line) 
which is scaled to remove the 11.3 /im PAH fea- 
ture. The spectrum on the top is the total ob- 
served spectrum while the dotted line shows the 
star-formation corrected spectrum (note the ab- 
sence of the 11.3 /Ltm PAH feature). 



order to explore any nuclear MIR emission in FR- 
I sources, we will here test this scenario by suc- 
cessively identifying contributions from different 
mechanisms to the MIR spectra. 

4.1. Contributions from star formation 

Due to the large slit width of IRS (3.7" and 
10.7" for SL and LL, respectively) contributions 
from the host galaxy (star formation and stellar 
emission) to the MIR continuum can be significant 
and even dominant. We need to correct our spec- 
tra for these host galaxy contributions in order to 
explore any additional nuclear MIR emission. 

Because emission related to star formation has 
intrinsically a red spectral slope it cannot be dis- 
tinguished easily from the synchrotron core emis- 
sion or AGN heated dust. Therefore, in order to 
securely identify any synchrotron or AGN dust 
components in the observed MIR spectra, emis- 
sion from star formation needs to be corrected for. 
As an indicator for star-forming contributions we 
here use the PAH features which are commonly 
observed in the spectra of star-forming galaxies 
and are a well known tracer of star-formation ac- 
tivity. Only if PAHs are detected do we correct the 
spectra for emission related to star formation. No 
significant PAH emission usually indicates negli- 



gible continuum contribution from star formation 
compared to other processes. 

In the case of detected PAH emission we 
chose to subtract appropriately scaled aver- 
age star-formation template spectra taken from 
.Smith et al., (2007.) . We used the template with 
the reddest MIR/FIR slope which was still in ac- 
cordance with the data after scaling, thus being 
conservative about any residual emission. The 
template was scaled to the 11.3 /im PAH feature 
which is strong, largely free of contaminations by 
atomic emission lines, and close to wavelengths 
where AGN related emission usually has a large 
contribution. The star-forming template was 
scaled in such a way that, after subtraction, the 
11.3 /im PAH feature is removed from our FR-I 
spectrum (see Fig. [5] for illustration). In Tables [2] 
and [3] we give estimates for the percentage contri- 
butions from star formation to the observed MIR 
fiux at 15 /im and 30 /im. 

While this strategy works well for many of our 
sources, subtraction of the star-formation tem- 
plate (scaled to the 11.3/im PAH feature) leaves a 
clear depression at the location of the 7.7 /im fea- 
ture for 3C83.1 and 3C272.1. This cannot be ex- 
plained by the shape of stellar emission and argues 
for a very low 7.7/im/11.3/im PAH ratio in these 
FR-I galaxies (Fig.[6|). Apart from the low PAH 
ratio both spectra can be explained by a combina- 
tion of processes in the host galaxy. We detect no 
significant nuclear continuum (but see §4.3). 

Such low PAH ratios have already been noted 
in local (ear l y-type) galaxies ( Kaneda et al. 2005t 
Smith et all 120071 : iKaneda et all l2008h . Sources 



with MIR continua dominated by star formation 
usually show no large variation in the 7.7 /im 
to 11 .3 /im ratio ( Brandl et al.ll200^ Smith et al 



120071 ) and similar star-formation dominated PAH 
ratios have been foun d in luminous ra dio-quiet 
and radio-loud AGN (jShi et all [2007bl ). In ab- 



sence of strong contributions from star forma- 
tion, spectra with unusual 7.7//m/11.3 /im val- 
ues can be observed. Interestingly, most sources 
with such lo w PAH ratios sho w low-luminosity 
AGN activity (ISmith et al.ll2007f ). Consistent with 
this, Sturm et al.l (|2006l ) have shown that IR- 
faint LINERs (LINERs without much star for- 
mation) show these unusual PAH ratios, while 
IR-luminous (Lir,/Lb > 1 with Lir being the 
8 - 1000 /im luminosity) LINERs have PAH ratios 
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Fig. 6. — The unusual PAH ratios in 3C83.1 {left) and 3C272.1 (right). The top panel shows the observed 
spectrum with the star-formation template scaled to the 11.3/im PAH feature. The bottom panel shows 
the results of the subtraction and the scaled spectrum of an early-type galaxy (NGC 1549). Approximate 
contributions from the synchrotron core (§4.4) are shown as long-dashed lines. 



similar to starburst galaxies. In the former cases 
it appears that the nuclear source powers most of 
the measured PAH emission. The harder spec- 
trum could potentially lead to the destruction of 
the carriers of the 7.7 /im band or to the preferen- 
tial excitation of the 11.3 /im PAH, thus explaining 
the o bserved relative decrease in the 7.7 iim fea- 



Kaneda et all l2008l and 



ture (|Smith et al.M2007 : 
references therein). This could in principle indi- 
cate that the PAHs in 3C83.1 and 3C272.1 are not 
excited by star formation but by (low-luminosity) 
AGN activity instead. Any corrections made to 
the spectrum by the subtraction of a star-forming 
template would thus remove MIR continuum emis- 
sion which is in fact powered by the AGN. 



4.2. Contributions from stellar emission 

Recall that we have already seen in the spec- 
tra (Fig.[T]) that contributions from stellar emis- 
sion of the host galaxy can be substantial, par- 
ticularly at shorter wavelengths. Most or all of 
our FR-I sou rces res ide in elli ptical galaxies (e.g. 
iGovoni et al 2000: Madrid et a l. 2006). The MIR 
spectra of early-type galaxies generally show blue 
col ors throug hout the IRS wav elength range (e.g. 
Br essan et al. 2006,; Kaneda et al. 2008) and even 
if moderate amounts of residual star formation 
are observed, the blue spectral slopes of the stel- 
lar emissi on can remain dom inant in the total 
spectrum ([Kaneda et al. I l2008h . This is supported 
by the fact that the SEDs of typical early-type 
galaxies without signs for considerable star for- 
mation usually show a blue IR slope that, due 
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to emission from cold dust, turns over into a red 
slope at w avclcnK ths beyond the IRS c o verage (e.g. 



Dale et al . 2005; Kaneda et alJ 120071 : iTemi et al. 
20071 ). Therefore, the stellar contributions in our 



MIR spectra will have opposite spectral slope com- 
pared to any underlying synchrotron core or warm 
dust emission. 

As shown by e.g. iBressan et"al] ([20061) the MIR 
spectra of non-active early-type galaxies can be 
well explained by including the dusty atmospheres 
and envelopes of AGB stars into stellar evolution 
models. The emission from dust produced in the 
atmospheres of AGB stars accounts for excess con- 
tinuum emission over the long-wavelength quasi 
Rayleigh-Jeans extrapolation of late-type stars as 
well as for the presence of silicate emission in the 
spectra of early-type galaxies. 

We estimate contributions from the stellar pop- 
ulation by using the observed spectrum of a quies- 
cent elliptical galaxy which is scaled to the blue 
part of the spectra (~5 — 8/im). For a tem- 
plate we chose the quiescent early-type galaxy 
NGC1549 which is fairly bright at MIR wave- 
lengths, does not show any PAH or atomic line 
emission, and has ful l IRS low-resolution coverage 
([Kaneda et al. l2008l) . We obtained the data for 



NGC 1549 from the Spitzer archive and extracted 
a spectrum fo l lowing the procedures described in 



Kaneda et all ([2008D . In Tables H and [3] we list 
estimated percentage contributions of the stellar 
component to the observed MIR flux at two wave- 
lengths. 

4.3. Spectral components in the MIR 

After the subtraction of the star-forming tem- 
plate in the way described above and after consid- 
ering the stellar host galaxy emission we recognize 
basically two types of spectra: 

(1) sources where the observed MIR emission 
can be accounted for by processes in the host 
galaxy. 

The spectra are either dominated by star forma- 
tion (3C31, 3C218, 3C264, 3C293) or the observed 
MIR emission can be explained by a combination 
of stellar and star-forming contributions (i.e. after 
the subtraction of the star-forming template the 
residual spectrum is largely consistent with stel- 
lar emission from an early type galaxy; 3C76.1, 
3C386; Fig.\T^. 



(2) objects which show significant residual emis- 
sion at wavelengths longer than ~ 8 /im with a red 
spectral slope. 

Because even star formation and stellar emis- 
sion combined cannot explain these residual con- 
tinua they must originate from other processes, 
for example warm (few hundred K) dust emis- 
sion or non-thermal synchrotron emission. In 
some of the sources where residual continuum 
emission at A > 8 fim could be identified we 
also detect PAHs (3C84, 3C120, 3C129, 3C270, 
IC4296, NGC 6251, E1821+643; Figs.[T3ll4|), but 
the equivalent widths of these features are gen- 
erally small. For many sources only red contin- 
uum emission at A > 8 /im (plus possibly stellar 
host galaxy emission) was detected (3C15, 3C29, 
3C66B, 3C189, 3C274, 3C317, 3C424, 3C465, BL 
Lac). No significant PAH emission can be seen in 
their spectra. 

As mentioned above in the case of 3C83.1 and 
3C272.1 we cannot without doubt attribute the 
PAH emission to star formation. We count these 
objects as host-dominated sources, an intepreta- 
tion which is supported by the shape of their MIR 
continuum (Fig. [HI), but we emphasize that we can- 
not be sure at this point if the PAH emission is ex- 
cited by star formation or (low-luminosity) AGN 
activity. 

4.4. Comparison with synchrotron predic- 
tion 

Since it has been proposed that non-thermal 
emission might be a major, if not dominant con- 
stituent of the MIR emission in FR-I sources, we 
also analyze the MIR spectra in the framework 
of their (nuclear) SEDs. As a comparison for the 
synchrotron core emission we chose to use well ob- 
served blazars and related objects because they 
represent nearly pure synchrotron emission. 

It has been shown that, globally over a large 
frequency interval, the strongly beamed core emis- 
sion of blazars can be well approximated by 
a parabolic functio n in logFi, vs. \ogv (e.g. 



Landau et al. 19861 ). Following these authors we 



here utilize a function of the form 



logF^ ^ C + [logiy ~ B)y2A 



(1) 



to represent the underlying beamed emission 
from the base of the jet (C is the log of the peak 
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Fig. 7. — Observed nuclear SEDs for sources in which significant residual continuum emission could be iden- 
tified. We used optical fluxes (Tab. 1), corrected for galactic absorption, nuclear NIR fluxes, and radio core 
flux measurements (arcsecond and sub-arcs econd scales). O nly optical/NIR and radio core measurements 
were used for the parabolic fits. Nuclear UV ( Chiaberge et al.| [2002) fluxes, also corrected for galactic absorp- 
tion, MIR/FIR, and sub-millimeter data (if available) as well as MIR spectra (this work) are overplotted. 
For sources with significant PAH features we have subtracted the star formation component (3C84, 3C270, 
IC4296, NGC6251; see text). 



flux, B is the log of the peak frequency of the 
parabola, and A represents a "curvature" param- 
eter). We note that a similar parabolic function 
can also be used to descri bed the SEDs of BL L acs 
in logv Fi, vs. logv (e.g. Nieppola et al.l 20061 ) . 

In order to construct the (nuclear) SEDs, pho- 
tometric data from the literature have been com- 
piled from optical through radio wavelengths. For 
the radio and optical we only focused on core emis- 
sion and we only compiled these data for sources 
with detected optical CCC components. In total 
our sampl e has 14 objects in com mon with the 



sample of Chiaberge et al. ( 19991 ). But for two 
sources no CCCs could be identified: for 3C293 



Chiaberge et al. ( 19991 ) could not estimate any 
CCC flux due to large scale dust structures in 
the host galaxy, and 3C424 showed a "radically 
different" nuclear behavior (i.e., the FWHM of 
the nuclear source is much larger than for other 
CCCs and it is spatially resolved). This leaves 12 
sources. In following works CCCs were also de- 
tected in 3C15 (R. Baldi as well as M. Chiaberge 
priva te communic ation), in 3C189 (ICapetti et al 



2002t ). in IC429 6 jBalmaverde et all l2006al) . md 



NGC6251 ()Chiaberge et al.l \200A Thus 



16 sources here have an optical CCC detected 
(Tab. 1). Compared to other CCC objects, 3C386 
stands out in the properties of its optical CCC 
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Fig. 7. — continued. 



which is exceptionally bright co mpared to the ra 



dio c ore. It has been suggested (|Chiaberge et al 



19991) . supported by the te ntative detection of a 
broad optical emission line ( Simpson et al. 19961 ) 



that this source shows additional optical flux from 
a Big Blue Bump component, i.e. showing the cen- 
tral core of a type-1 AGN. However, it seems more 
likely that the foregro und star which falls right 
on top of the nucleus (iLvndd 1971 : Madrid et al 



20061 : Buttiglione et all 2009 ) mimics an optical 
CCC because we do not see any signs for such a 
type-1 AGN in the MIR. In fact, after subtracting 
some emission due to star formation the residual 
spectrum can be explained well by a quiescent 
early-type galaxy (Fig. [13]). In addition, new op- 
tical spectroscopy does not confirm a broad Ha 
emission-line component ( Buttiglione et al.]l2009f l . 
Because we cannot tell whether or not 3C386 it- 
self shows a compact optical core we here exclude 
this object from further analysis. This leaves 15 



sources where optical CCCs are securely detected. 
Eleven of the 15 objects also show residual MIR 
emissio while the spectra of four sources are 
dominated by processes from the host galaxy (stel- 
lar emission, star formation). 

We used the parabolic function to fit only 
the radio core data and the optical CCCs which 
are both claimed to be synchrotron emission 
from the same source re.g.lChiaberge et al.|[l999t 



Hardcastle fc Worralll [2000) 7" For these data we 
know that they are nuclear and fairly robust while 
other measurements might still be affected by ei- 
ther intrinsic absorption (UV) or aperture effects 
(MIR, FIR, sub-mmB 



^Recall that if our methodology is correct the presence of 
residual emission means they contain AGN heated dust 
and/or prominent synchrotron contributions, i.e. a nuclear 
source of MIR continuum emission. 

^Because the core flux data were compiled from the litera- 
ture they are generally not obtained simultaneously which 
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4-.4-1- Sources with residual MIR 



emission 



In the following we will discuss the SEDs on an 
object to object basis. References to flux measure- 
ments used for the SED plots (Fig. [7]) in addition 
to those already given in Tab. [T] are also noted. In 
Tables [2] and [3] we give the percentage contribu- 
tions of the estimated non-thermal component to 
the observed MIR flux. 

3C15 This source shows MIR continuum 
emission which is not due to either star forma- 
tion or stellar emission (Fig.[3|). Using optical and 
NIR CCC measurements (provided by R. Baldi as 
well as M. Chiaberge, private communication) and 
radio core data the non-thermal core component 
clearly underestimates the observed MIR flux. 



can result in variability becoming important for the quality 
of the fits. However, in most cases all the multi-epoch data 
can be fitted well and any scatter present due to intrinsic 
variability seems relatively small. 



However, we poin t out that 3C15 ha s a strong 
arcsecond-scale jet ( Leahv et alj 1997) which is 



also detected in the optical 
and in X-rays (jKataoka et al 



Martel et all 119981 ) 
200^. Since the 



extended jet is included in both IRS slits (SL slit 
oriented along the jet axis, LL slit oriented per- 
pendicular to the jet axis) it can potentially make 
additional non-thermal contributions to the MIR 
spectrum. The inspection of the 2-D spectral im- 
ages from IRS does not reveal any secondary com- 
ponent in the SL profile but the spectrum appears 
dominated by a single source (the size of the jet is 
almost comparable to the resolution of Spitzer at 
IRS wavelengths). 

At 8.4 GHz the total radio flux of the arcsec- 
ond jet is greater than the core flux by a facto r 
of ~10 (|Leahv et al.lll997tlHardcastle et al.lll998f) . 
Assuming a typical jet spectral index of a ~ —0.6 
(jBridle fc Perlevlll984j ) this would lead to a non- 
thermal jet contribution of 2 mjy and 3 mjy at 
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Fig. 8. — Same as in Fig. [7] but for sources which are dominated by processes in the host galaxy (star 
formation and/or steUar emission). Here, we show the observed MIR spectra. 



15 ^m and 30 fj,m, respectively, in the case of a sin- 
gle power law. For this conservative estimate we 
used a spectral index which is quite flat for a jet. 
In fact, the radio-optical spectral inde x of the X- 
ray b rightest jet knot ("Knot C", Kataoka et al 



20031) is argued to be steeper than our chos en value 
(|Kataoka et"ani2003l : [Pulwich et al.ll2007^ . 



Thus, even under the assumption of a conser- 
vatively flat spectral index for the jet the total 
non-thermal contribution (core -I- jet) to the MIR 
emission within the IRS apertures is still a factor 
of > 2 smaller than the total observed emission. 
In addition, the spectrum itself shows a change 
in slope around ~ 20 /im (which could be identi- 
fied with a 18 /im silicate feature) which argues 
for thermal contributions to the total flux. 

Considering the appearance of the spectrum 
(including the potential 18 /im silicate feature), 
the absence of notable PAH emission, and the dif- 



ference between total flux and (conservatively es- 
timated) non-thermal contributions we therefore 
argue that 3C15 shows thermal MIR emission from 
warm nuclear dust, possibly heated by an AGN. 

We also note that iRinn et al.l (|2005[ ) ob- 
serve the X-ray core of 3C15 to be fairly weak 
(L2-iokeV ~5 X 10''°erg/s intrinsic) but absorbed 
by a significant column density of A '^h ~ 9 x 
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22 , 



Re ferenc e s: iDicken et al. (12008 ) ; 



iHardcastle et al.l (|l998l ): iMorganti et al.l (|l997t ) 

3C29 The parabolic fit matches the radio 
data quite well but greatly overpredicts the nu- 
clear UV measurement. Assuming that this is due 
to extinction, dereddening of the UV flux in or- 
der to match it with the synchrotron fit would 
also result in increasing the optical flux somewhat. 
Such a synchrotron SED could in principle ac- 
count f or the small excess in t he MIR spectrum. 
In fact, IChiaberge et al.l (|2002f ) note that the UV 
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Fig. 9. — The star-formation corrected spectrum 
of 3C270 is fitted with the spectrum of the el- 
Hptical galaxy NGC 1549 (dotted fine), the syn- 
chrotron component from the core (dashed Une), 
and two black bodies (dot-dashed line; see text for 
details). 



point in 3C29 might be absorbed by a thin, ex- 
tended dust lane. However, the mismatch of the 
UV measurement could also arise from variability. 
While the optical and UV data were obtained more 
than 5 years apart, the difference between the 
synchrotron estimate and the actual UV measure- 
ment corresponds to a factor of ~ 30 in flux. This 
seems high to be solely due to variability. Con- 
sidering that some absorption in the UV/optical 
bands is present, that the slope of the MIR spec- 
trum agrees reasonably well with the synchrotron 
fit, and that no host corrections were applied the 
MIR emission in this sour ce is likely to be mostl y 
non-thermal. Ref erences: iMorganti et"al] (|l997l ): 



Ricci et al. 



tOO&) 

3C66B Here all data points (including the 
UV and MIR) agree very well with the parabolic 
fit. Thus, there does not seem to be sig- 
nificant extinction present in this object and 
the residual MIR emission is very likely to be 
non-thermal. At the shortest wavelengths of 
the MIR spectrum the stellar population of 
the host galaxy c an be seen. 3C66B has an 



arcsecond-scale jet (jHardcastle et al.lll996l) which 
is also seen in the op t ical (iButcher et al. 19801 : 



Macchetto et al.l 1991 : Perlman et al 



X-rays (jHardcastle et al.ll200l[ ). and even at MIR 
wavelengths a detection with ISO is claimed 



(jTanslev et al. l2000h . However, unlike the case 
of 3C15, the jet i n 3C66 B is weaker than the core. 



Hardcastle et al.l (|2001l ) fit the radio through X- 



ray jet spectrum with a broken power-law model 
where the steepening from a ~ —0.5 to a —1.35 
occurs in the infrared. The radio-IR spectral 
slope agrees reasonably with the 14.5 /xm MIR 
j et flu x of ~1.7mJy estimated bv iTanslev et al 
( 2000l ). Considering these values the jet contri- 



bution in the MIR is only < 50% of that of the 
core. Thus, the inclusion of the jet in the IRS 
spectral aperture does not alter our conclusion 
for the dominance of non-thermal MIR emission 
in nucleus of 3C66B. Refe rences: Ijackson et al 



Kharb et all ^QQM) 



(Il993l): IXu et all (l2000l): lOuillen et all (l2003h : 



3C84 While this source is strongly core dom- 
inated in the radio, a clear dust emission bump is 
observed which reaches from the FIR well into the 
MIR. The SED is very complex and suffers from 
resolution effects (sub-mm vs. VLA vs. VLB A), 
which reduce the quality of the synchrotron core 
fit. In addition the VLB A data suggest self ab- 
sorption in the core making the SED asymmetric 
which cannot be fitted well with a parabolic func- 
tion. Therefore, the fit was limited to the optical 
core and the radio core measurements excluding 
the 5 GHz data where the downturn in the SED 
is already quite prominent. Besides the resolution 
effects and the self absorbed core, variability might 
also be important and the fit cannot be more than 
a rough approximation of the synchrotron core 
component in this case. However, the thermal 
MIR and FIR bump is very pronounced in 3C84 
and variations in the underlying synchrotron com- 
ponent are almost negligible at these wavelengths. 

The MIR spectrum itself does not show sig- 
nificant contribution from stellar emission, only 
minor indications for star formation, and overall 
has a spectral appearance quite similar to "ther- 
mal" spectra of AGN powered sources. The fact 
that the spectrum is on average almost an order 
of magnitude brighter than the underlying non- 
thermal core emission strengthens the argument 
for thermal dust emission i n 3C84, powered by a 
central engine. In addition. iHoet alJ (Il997bh find 
broad emission-line components in Ha and other 
ermi tt ed lines. Re f erenc es 
1991: iQuillen et al l (l2003h 



Knapp fc Patten 



Haas et all (|2004l ): 



Stickel et all (|2n04f l: TTavlor et al.l (|2006[ ) 
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3C189 In the observed MIR spectrum this 
source is very similar to 3C66B (Fig.[T]). Also 
in the context of the SED we see that the fit 
to the synchrotron core matches with the MIR 
emission in shape as well as in flux. An IRAS 
100 /j,m flux measurement indicates the presence 
of cool (possibly e xtend ed ) dust. Re fe rences : 
Bridle & Fomalonll ("lEITSh: Ixnapp et all (|l990l ): 
Giovannini ct al. (199^) 

3C270 Including nuclear NIR measurements 
as well as (large aperture) FIR and (sub)-mm 
data, the SED has very good coverage. We see 
that the radio core data, the NIR, and the (sub)- 
mm data fit well with the synchrotron predic- 
tion, arguing for their common non-thermal ori- 
gin. The MIR spectrum, however, shows a very 
clear excess over the underlying non-thermal emis- 
sion. Not only does the excess in MIR flux indicate 
additional thermal emission from warm dust in 
this source, but also the shape of the Spitzer spec- 
trum makes its nature as non-thermal emission 
very unlikely. The pronounced change of slope 
around A ~ 20 /im rest wavelength (Fig.[5|) strongly 
suggest thermal emission from warm dust to be 
present in this source. 

In fact, after subtracting the star-forming tem- 
plate from the spectrum the MIR continuum can 
be well fitted as a combination of the underlying 
synchrotron contribution and the elliptical tem- 
plate for the host galaxy, plus a black body with 
a temperature of ~ 200 K (Fig.[Hlll|- This emission 
from warm dust is not powered by star formation 
because this would result in accompanying PAHs. 
At 15 /im the continuum emission of this resid- 
ual thermal MIR component has a luminosity of 
i^ii^,i5Mm ~4 X 10"'-^ ergs~^ (see also Tab. [3]). 

In the optical, iBarth et al. report a 

broad Ha line in polarized light. This (tentative) 
detection further supports the presence of a classi- 
cal hidden AGN. Moreover, a geometrically thick 



nuclear dust disc has been observed in silhouette 
(jjaffe et al.l Il993l ). We also note that the opti- 



'^The fit to the star— formation corrected spectrum has only 
3 variable parameters: the scaling of the elliptical template 
as well as the temperature of the two black bodies. Be- 
cause the synchrotron component was determined indepen- 
dently from radio/optical core measurements its flux was 
held fixed. We note that the colder black body is not well 
constrained but including it improves the fit for A > 20 fim. 
However, the presence or absence of this second black body 
leaves the temperature of the warmer component virtually 
unchanged as it is strongly constrained by the continuum 
between 10 and 20 ^m. 



cal emission-lines in the nu clear region have flux 
ratios expected for Seyferts (jFerrarese et al. I ll996l: 
see §4.8). 

Several studies are in broad agreement on the 
basic X-ray properties like the moderate intrinsic 
luminosity of L2-iokcV~l x 10^^ erg/s and the 
moderate a bsorption on the orde r of A^h 5 x 
10^^ c m-^ dChiaberge et all I2OO3I: iGliozzi et al 



2003 



Sambrun a et all 120031; iRinn et all l2005t 

Zezas et al. 200?C The origin of the observed 
X-ray emission is, however, co n trover sially dis- 
cussed: While IChiaberge et al.l (|20031 ) interpret 
the X-ray power-law c omponent as due to the 
jet, Gliozzi et al. ( 2003f ) argue in favor of an ac- 
cretion flow as the most likely source for the 
bulk of the X-ray emission. On the other hand, 
Zezas et al.l (|2005l ) state that both processes are 
able to explain the observed X-ray properties. 
A Fe Ka line is (marginally ) detected with an 
equivalent width of ^ 230 eV dGliozzi et al. 20031: 
Sambruna et al. I I2OO3I: iRinn et al. I I2OO5I) which is 
cons istent with the upper limits from other stu dies 
(e.g. IChiaberge et al.ll2003t IZezas et al.ll2005l) . 



While the warm dust might represent the same 
nuclear absorbing material detected in X-rays 
we note that due to the large Spitzer slit width 
our spectra include the prominent ^ 300 pc nu- 
clear dust disc. In principle the 200 K dust 
observed in the MIR could correspond to this 
dust disc. Because in situ star formation as 
a heating mechanism can be largely excluded 
(no PAHs) another source of sufficient energy 
needs to be present in the nucleus of 3C270. 
Radiation from a (possibly hidden) active nu- 
cleus is an obvious possibility. However, in the 
case of heating by the low-to-medium luminos- 
ity AGN in 3C270 simple energy budget argu- 
ments strongly suggest the ^ 200 K dust to re- 
side on much smaller scales than ~ 300 pc. Ref- 



erences: Jones & Wehrle 


( 

); 


1997t): Caoetti et al. 


(I2OOO): Jones et al. (2000 


Quillen et all (2003): 


Haas et al.l (2004): Kharb et all (2005) 



3C274 Interestingly, the slope of the Spitzer 
spectrum appears flatter than what one would ex- 
pect according to the synchrotron core fit. After 
accounting for host galaxy light using the ellip- 
tical template the MIR continuum alone can be 
well fitted with a single power law of a '^—0.80 



19 



for wavelengths shorter than ~ 25 /im. At longer 
wavelengths the spectral slope steepens with re- 
spect to the power law which mi ght be due to a n 
additional dust co mponent dPerlman et al. 20071) . 
As discovered by IWhvson g fc A ntonuccil ( 2001 ) 
and later discussed by iPerlman et alj ( 2001 ) 
and Whysong fc Antonuccin 20041 ) . an unresolved 
~ 0.5 " MIR core is present in 3C 274. The MIR 



spectr um is discussed i n detail bv iPerlnian et al . 
(I2OO7I ). Refer ences: iKnapp fc PattenI h99W. 



IXu et al .l(l2000l):lNagar et al.l(l200l[) :'Quillen etall 
(2003); Haas et al.l (l2004l);lListerfc H oman (200£ 



Pcrlman et al 



Whvsong fc Antonuccil ()2004D 



( 20011. I2OO7I): IShi et all (|2007a 



IC4296 For IC4296 the synchrotron fit falls 
very close to the MIR spectrum, but only for a 
very narrow wavelength interval (A ~ 20 — SO/im). 
In this source a star-forming template has been 
subtracted because PAHs are clearly detected and 
the spectrum seems to be dominated by stel- 
lar emission shortwards of A~10/im (Fig. [14]). 
Above 25 /im a considerable steepening is ob- 
served in the spectrum suggesting the presence 
of a thermal bump in the SED due to cooler dust 
which is supported by FIR measurements. De- 
spite the influence of processes in the host galaxy 
at some IRS wavelengths, this source has a very 
strong non-thermal contribution of greater than 
> 60% at 15 /im and 30 fim (Tab.[2]). Ahhough the 
synchrotron signature is weaker than for sources 
like 3C66B or 3C189 it manifest s as clear resid- 
ual MIR emission. Interestingly, IPelleg rini et al 



([2003) argue that the X-ray core in IC4296 
(L2-iokcV 1-6 X lO'^^erg/s intrinsic) originates 
from the jet and the X-rays are observed to be 
only m ildly absorbed (A^h ~ 1x10^^ cm~^) Refer - 
ences: iKilleen et al.l Il986l): iKnapp et alJ (Il989l): 



the reasonable match in flux (combined with the 
absence of star-formation tracers) strongly sug- 
gest the dominance of non-thermal emission in 
the MIR. Spitzer FIR measurements at 70 /xm sug- 
gest the presence of cold d u st in this obje c t. Ref - 



Zhao et al.l (Il993l) : IVenturi et al.l (|200ol ): 

Quillen et all (|2003L 120081 ) 



NGC 62 51 Using UV, optical and NIR mea 



surements (jChiaberge et al.l |200^ and multr 



freque ncy radio core measurements (VLBI. lEvans et al 



2005br ). a well sampled nuclear SED can be con- 



structed for NGC 6251. As for other sources, vari- 
ability might be an issue for the non-thermal core 
measurements and in fact variability with a fac- 
tor of ~ 2 over a few years has been reported for 



the nu clear (VLBI) radio component (| Evans et al 
2005bl ). However, the radio measurements shown 



iMorganti et al.' ('l997l):IVe nturi et al.l (|2000l ) dPellegrini et 

(^.2003); Ricci et al. (2006) 

3C317 For this source the MIR spectrum 
falls very close to the synchrotron flt in flux. 
The synchrotron curve overestimates the nuclear 
UV measurement somewhat and a similar situa- 
tion of absorption intrinsic to the source or vari- 
ability as in 3C29 co uld be imagined. In fact, 
Chiaberge et al.l ([2002) reported a factor of ^ 10 in 



here were obtained simultaneously and the fit also 
agrees well with most of the multi-epoch HST 
measurements. For this source the MIR spectrum 
shows a clear excess over the synchrotron core fit, 
even after the correction for some minor star for- 
mation as traced by PAH emission of small equiv- 
alent width (Fig.[T]). The small contribution from 
star formation, the absence of signs for significant 
contributions from stellar emission, and the over- 
all shape of the spectrum (including the strong 
silicate emission features) strongly suggests a con- 
siderable component of thermal MIR emission in 
the nuclear regions of NGC 6251. The spectrum 
also looks quite different from other sources which 
are dominated by non-thermal emission in the 
MIR (e.g. 3C66B, 3C189 in Yig^. 

At 30 /im, where the strong silicate emission 
features will have only a minor impact on the 
continuum luminosity (e.g. lSpoon et al.ll2007l ). we 

_Li|easure uL„ = 5.6 x lO'*^ ergs~^ in the spectrum. 

S^tl the other hand, the luminosity of the syn- 
chrotron component at this wavelengths is about 
vLi,^2.2 X lO'*^ ergs^^, which is consistent with 
the vL^^2.b X 10*^ ergs^^ gi ven by the syn- 



chrotr on self-Compton SED fit of lChiaberge et al 



variability for the nuclear UV point source. Thus, 
the discrepancy of the CCC measurements (and 
the offset in MIR flux) is likely to be caused by 
variability. The similarity in spectral shape and 



()2003l) . We note however, that the factor of ~2.5 
luminosity difference between the synchrotron 
core and the MIR spectrum is of the order of 
the measured radio core variability. In principle, 
the excess of the spectrum could be accounted for 
by variability of the underlying synchrotron core 
component. Although this would change the ther- 
mal to non-thermal ratio in the MIR (with the 
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Table 2: Estimated contributions to the MIR spectra for sources without detected AGN dust components. 



Object 


Fi5 (rest frame) 


F30 Mm (rest frame) 






total^ 


SF 


stars 


synch 


total^ 


SF 


stars 


synch 


comment 




mJy 


% 


% 


% 


mJy 


% 


% 


% 




sources without optical core measurements 


3C76.1 


1.72 


21 







3.53 


30 







low S/N spectrum 


3C129 


3.26 


17 


46 




11.44 


14 


6 






3C218 


4.49 


100 







15.91 


100 







SF dominated spectrum 


3C293 


20.10 


100 







57.98 


100 







SF dominated spectrum 


3C386 


1.97 


16 


74 




1.78 


52 


37 




foreground star mimics CCC 


3C403.1 


< 1.34 








< 1.53 








not detected 


3C424 


1.52 




12 




5.06 




2 






sources with optical core measurements 


3C29 


2.48 




15 


58 


8.16 




2 


42 


peak-up spill-over^ , UV abs. 


3C31 


13.99 


87 





13 


29.30 


86 





14 


SF dominated spectrum 


3C66B 


4.76 




21 


79 


8.72 




5 


85 




3C83.1 


4.08 


18 


64 


12 


5.37 


40 


22 


22 


low PAH ratio 


3C189 


4.80 




20 


71 


8.73 




5 


79 


cold dust 


3C264 


9.11 


52 





48 


23.88 


67 





33 


SF dominated spectrum 


3C272.1 


28.08 


30 


57 


13 


46.70 


53 


15 


15 


cold dust, low PAH ratio 


3C274 


43.25 




7 


93 


81.38 




2 


98 




IC 4296 


10.02 


6 


30 


69 


19.88 


9 


7 


68 


cold dust 


3C317 


2.90 




16 


50 


5.86 




3 


57 




3C465 


3.08 




19 


68 


7.47 




3 


59 


cold dust 


BL Lac 


246.45 







100 


381.92 







100 





"Flux measured from the observed MIR spectrum. 

''Flux of the synchrotron component shown in Figs. [7l and |81 

'^Due to the peak-up spill-over the spectral slope and flux at short wavelengths might be compromised. This adds some 
uncertainty to the scaling of the stellar template. 

Note. — The contributions of the individual components do not necessarily add up to 100%. This is due to the nature 
of our appraoch which uses fixed templates in a step-by-step manner in order to identify and correct successively for their 
contributions to the observed spectra. This naturally leaves room for variations between individual objects which can manifest 
as residuals or slight overcorrections. Also, components not considered might have an effect (e.g. cold dust, particularly at 
30/xm). In addition, for sources with significant contributions from non-thermal emission variations in the nuclear components 
and obscuration for UV/optical data might result in some mismatch and residual emission. See text for comments on individual 
sources where this might be the case. 



non-thermal part potentially even dominating at 
times) the appearance of the MIR spectrum still 
argues for the existence of warm nuclear dust in 
NGC6251. 

Ferrarese fc Ford ( 1999( ) report the tenta- 



tive detection of a broad component in Ha 
(this result seems to be strongly dependent 
on the assumed modelling of the narrow lines 
though). In X-rays the power-law component 
in NGC6251 shows only very little absorption 
(A^H^S X 10^"cm~^) and has an intrinsic lumi- 



nosity ofL^ 



^4.8x lO''^ erg/s (|Gliozzi et al 



2004 lEvans et al.ll2005bh . The additional detec- 
Ka line (EW ~220eV and 



Fe 



tion of a broad 

FWHM~0.6keV) led iGliozzi et al.l (|2004D to ar- 



gue in favor of the presence of a standard accretion 
disk. In their study the base of the jet is ruled out 
as the sole origin of the nuclear X-ray emission. 
Evans et al. ( 2005lJ ) on the other hand question 



the significance of the Ka detection and favor the 
jet as the source for the bulk of the observed X- 
ray core emission. However, they do not exclude 
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some co ntributions from an acc r etion flow. Ref- 
erences: KJolombek et al. ' ('1988'); ChiabcrEC et alj 
(120031); [Quillen ct al 



Evans ct al. (2GG5b|) 
spectrum falls very 



(2003) 

3C465 The observed 
close to the synchrotron core fit. Assuming that 
the observed nuclear UV component is also non- 
thermal, significant absorption seems to be present 
in this source. Alternatively, variability of a fac- 
tor of ~ 10 could cause this UV flux difference. As 
noted before, both processes can in principle ac- 
count for the remaining offset between the MIR 
spectrum and the synchrotron fit. A steepen- 
ing of the spectrum is observed for wavelengths 
larger than ^ 30 /im indicating a thermal bump at 
FIR wavelengths. This is suppo rted by IRAS 60 
and 1 jum fluxes. Refere n ces: iGolombek et al 



'Venturi et all (|l995l) ; lQuillen et al.l (|2003l ); 

Jetha et a l. (20061) 

BL Lac U sing quasi-simultane ous flux mea- 
surements from lLandau et al I (|l983l) we construct 
an SED of this object, which is here used as an ex- 
ample for a strongly aligned, highly beamed FR-I 
source. Note that although no HST data were used 
the (ground-based) optical fluxes will be domi- 
nated by the core emission. The measurements 
can be well fit by the parabolic function and even 
though the MIR spectrum was obtained at a differ- 
ent epoch, it fits extremely well with the estimate 
for the synchrotron core emission. The spectrum 
can be fitted with a power law with a ~ —0.7. 

( 1983h 



References: iLandau et al 



4-4-^- Sources without residual MIR emission 

In the following we will present the sources with 
detected CCCs but where no significant residual 
MIR continuum emission can be identified. The 
observed spectra can be explained by star forma- 
tion and/or stellar emission. The SED plots are 
shown in Fig. [51 

3C31 This source is strongly dominated by 
star formation in the host galaxy (Fig.[T]). Its 
observed spectrum is almost an order of magni- 
tude brighter than then underlying non-thermal 
component. Interestingly, the SL aperture for 
the MIR spectrum is roughly consistent with the 
size of the mol ecular disk observe d in the center 
of this source (|Okuda et al.ll2005^ . IRAS 60 /im 
and 100 /tm measurements indicate subs tantial 
amounts of colder FIR dust. References: iBurch 



(|l977l ); IXu eTall (|2000l ); lEvans et al.l (|2005al) 

3C83.1 The few existing data points for 
the core emission are well approximated by the 
parabolic fit. The total light spectrum is domi- 
nated by processes in the host galaxy and lies well 
above the synchrotron fit. The subtraction of the 
star-forming template leaves a negative residual 
at the location of the 7.7 PAH feature (Fig. [6]). 
As discussed above, this cannot be accounted for 
by the stellar spectrum of the host galaxy and ar- 
gues for a low 7.7 /im to 11.3 /im PAH ratio which 
might in fact be powered by low-luminosity AGN 
activity instead of star f o rmati o n (see ji4. 1 ). Re f- 



O'Dea fc Owml (|l986h ; IXu et"aD (|l999l ) 



3C264 While the observed spectrum of this 
source is dominated by star formation, it falls 
much closer to the synchrotron component than 
in 3C31. From Fig.[8]we see that the non-thermal 
contribution might not be negligible (see also 
Tab. (2). However, without prior knowledge of this 
fact the source spectrum can be fitted by star for- 
mation only due to the wide variety in the MIR 
spectral propert i es of star-forming galaxies (e.g. 
Dale et al.l [2OO5I ; iBrandl et all l2006t Ismith et al ' 



20071) . 



We note th at this source has an arcsecond-scale 
radio jet (e.g. Lara et al. 20041) which is also de- 



tected at optical wavelengths ( Crane et al. 19931 



Baum et al.lll997 ). However, in the radio the core 



is stronger than the total jet emission and ha s 
a significantly flatter spectrum (jLara et al.ll2004r ). 
While the jet emission steepens even more at op- 
tical wavelengths the contributions in the MIR 
are negligible c or npared to the cor e . References: 



Xu et al.r(|2000( ); iLara et all (|2004[ ); iKharb et al 



(|2005D 

3C272.1 The unusual sharp turn from a 
slightly blue slope into a steep, red slope observed 
in the MIR spectrum of this source is supported in 
the framework of the SED where several FIR mea- 
surements indicate the presence of a strong ther- 
mal bump from cool dust. The spectrum lies well 
above the core synchrotron fit at almost all wave- 
lengths and comes closest near A ~ 25 /im where 
the slope change takes place (Fig.[6|). The com- 
bination of strong PAH emission and a very low 
PAH 7.7 /im to 11.3 /im ratio results in a promi- 
nent depression at the location of the 7.7/im fea- 
ture after subtraction of the star forming tem- 
plate. As discussed above (§4.1) the unusual 
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Table 3: Estimated contributions to the MIR spectra for sources with detected AGN dust components (see 
notes to Tab. ID). 



Object 


F 


L5 f-LUl 


(rest frame) 




F 


30 fim 


(rest frame) 








total 


SF 


stars 


synch 


AGN^ 


total 


SF 


stars 


synch 


AGN^ 


comment 




mJy 


% 


% 


% 


% 


mJy 


% 


% 


% 


% 




3C15 


5.52 




19 


14 


67 


9.15 




5 


18 


77 


peak-up spill-over'' 


3C84 


1307.88 


1 





12 


86 


3859.70 


1 





9 


90 




3C120 


345.79 


2 







98 


653.08 


11 







89 


type-1 AGN 


3C270 


27.45 


7 


21 


11 


61 


52.14 


27 


5 


16 


52 




NGC6251 


27.88 


3 





29 


68 


49.68 


5 





33 


62 




E1821+643 


333.78 


2 







98 


540.74 


9 







91 


type-1 AGN 



"The contribution of the AGN component (which dominates in these objects) was determined as the total flux minus the other 
three contributions considered here. 
'See comment on 3C29 in Tab.0 



PAH ratio might point toward AGN excitation 
instead of star formation but we cou nt this source 
here as host dominated. Ref erences: Bower et al 



Haas et al.l ( 2004 ) : Leeuw et al. I (l2004l) : iDoi et al 
(l2005l) 



(l2000l): [Leeuw et a .1 toodi): iNagar et all (l200lh: 



4.5. Sources without CCC measurements 

In §4.4 we study the contributions of the syn- 
chrotron cores to the MIR emission in our sources. 
However, this test is limited to objects where opti- 
cal CCC measurements are available. When using 
the CCG fluxes in combination with radio core 
data, a nuclear SED can be constructed which 
bridges the MIR wavelengths. This way a fairly 
robust estimate of the non-thermal component 
can be achieved which helps to distinguish be- 
tween thermal and non-thermal contributions to 
any residual MIR emission (i.e. emission not re- 
sulting from star formation or stellar processes). 

Where no CCC measurements are available, a 
reliable estimate of the synchrotron core contribu- 
tion to the MIR emission is not possible. How- 
ever, the overall shape of the spectrum might still 
provide indications on which process - thermal or 
non-thermal - dominates the MIR emission. Out 
of the 25 sources considered in this paper, 7 ob- 
jects do not have CCC measurements available. 
We here briefly comment on these seven sources: 

3C76.1 The MIR spectrum of this source can 
be well explained by a combination of star forma- 



tion and stellar emission. No significant residual 
continuum can be identified. However, this spec- 
trum suffers from low S/N which limits the relia- 
bility of any conclusions regarding the MIR con- 
tinuum constituents. 

3C129 After correcting for some contribu- 
tions from star formation and accounting for stel- 
lar emission from the host galaxy, 3C129 shows 
strong, red continuum emission longwards of 
^15^m (Fig.fTH]). Inspecting Fig. [T2] reveals that 
we could in principle have used a star-formation 
template with a steeper MIR/FIR slope which 
would reduce the level of residual continuum. 
However, given the lack of any strong optical emis- 
sion lines, the general appearance of the optical 
spectrum, and the clear domina nce of an old stel- 
lar population (as determined bv lButtiglione et al 



.2009 ). it seems unlikely to us that this source 
harbors the stronger star formation which could 
power the redder MIR/FIR colors of a different 
template. However, while residual nuclear con- 
timmm emission is observed in 3C129 we cannot 
discriminate between a thermal or non-thermal 
origin because we lack optical CCC measurements. 
But we note that the residual MIR spectrum shows 
some concave upward curvature (Fig.ll3p and is 
significantly different in appearance compared to 
objects like 3C66B or 3C189 in which the MIR 
continuum is dominated by non-thermal emission 
(Figs. [TT7| . This argues in support of the thermal 
origin of the residual emission. 

3C218 The MIR spectrum of Hydra A can 
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Fig. 10. — Curvature parameter plotted against the logarithm of the peak frequency for the core SED 
fits in (left) and lyFjy {right). Open symbols represent FR-I sources, filled symbols the blazars from 
Landau et al.l The solid and dashed lines show linear least-square fits to the blazar and FR-I data 



points, respectively. Note that a smaller absolute number for the curvature parameter corresponds to a more 
strongly curved SED. We marked the sources BL Lac, 3C84, and 3C274 which are quite different from the 
lobe-dominated objects of our core sample but show extended radio emission of FR-I type. 



be well explained by star formation only, per- 
haps with some minor contributions from stel- 
lar emission at the shortest MIR wavelengths 
(Fig. [2]). Considering the equivalent width of the 
PAH features the continuum at A > 20 /im is in 
agreement with what can be expected for nor- 
mal star forming galaxies. The presence of a very 
young s tellar population is also seen in optical 
spectra iM elnick et al.lll997tlAretxaga et allfioOl : 
Wills et al.,.2004) . 

3C293 The MIR spectrum of 3C293 is 
strongly dominated by star formation and no 
residual nuclear continuum can be identified 
(Fig.[T]). In these properties it is quite similar 
to 3C31. 

3C424 For 3C424 no significant contribu- 
tions from star formation can be identified in the 
MIR spectrum and a weak and noisy, but securely 
detected, continuum can be seen. 

3C120 and E1821+640 are broad-line AGN, 
which have MIR spectra typical of their class, dis- 
playing strong silicate emission features and high 
ionization emission lines (Fig. [2]). 

4.6. Comparison v^rith blazars and related 
objects 

It has been suggested that most (but not 
all) FR-I radio sources are the misaligned coun- 



terparts of BL Lac objects (IKollgaard et al 
19921: lUrry &: Padovanilll995HCassaro et allllQea 
Antonuccill2002a[ ). As outlined in §4.4, the SEDs 



of these strongly beamed, highly core dominated 
sources are well described by the same parabolic 
function we utilized to fit the cores of the FR-I ob- 
jects (which is here used just as a fitting function 
for descriptive purposes). 

In order to allow the comparison between the 
blazars and the FR-I cores we used the q uasi- 
simultaneous data from lLandau et al" ( 198d ) to fit 
the blazars in the same way as the FR-Is. How- 
ever, we limited the parabolic fit for the blazars 
to the radio (2 cm through 20 cm) plus optical R- 
band measurements, thus maximizing the compa- 
rability with the FR-I fits (the wavelengths of the 
nuclear FR-I data used to constrain their fits were 
very similar). Overplotting all the blazar data 
onto the fits obtained using only the limited data 
points shows that this approach yields a reason- 
able representation of the curved SEDs. 

We will now compare the resulting fit parame- 
ters for peak frequency and SED curvature of the 
(putatively) misaligned sources with the respec- 
tive data of their highly aligned (supposed) sib- 
lings. Moreover, we performed the SED fits not 
only in F^ but also in vF^, using the same parabolic 
function as given in equation (1) and using the 



24 



same set of nuclear data. 

As already shown in e.g. Landau et al. ( 19861 ) 



Table 4: 24 /^m and [O III] emission-line fluxes. 



a correlation of peak frequency and curvature in 
exists for blazars. For FR-I sources a very 
similar behavior is detected (Fig. [TO|) . However, 
while both types of objects cover about the same 
range in F,^ peak frequencies there is a system- 
atic offset for the FR-Is towards more strongly 
curved SEDs. A comparable result has been ob- 
tained using broad band spectral indices of FR- 
Is and low-energy peaked BL Lacs. These two 
types of objects overlap in their radio/optical 
power-law slopes but there is a clear tendency 
for FR-Is to have steeper a^o tha n the BL Lacs 
(iHardcastle fc Worrall.2000; Trussoni et al.|[2005 
Balmaverde et al.ll2006al ). Because for a parabola 
peaking at radio wavelengths (which is typically 
the case for the sources considered here) a steeper 
radio/optical single power-law slope corresponds 
to a higher curvature, our results shown in Fig. 1101 
are consistent with these previous findings. 

In the blazars also show a correlation be- 
tween peak frequency and curvature, but now the 
FR-Is seem to continue this correlation at lower 
peak frequencies. However, we see no trend in 
the blazar distribution with FR morpholog 
particular, blazars with extended radio luminosi- 
ties typical for FR-I objects do not necessarily 
concentrate towards the "classical" FR-I objects 
in Fig. llOl but seem to be randomly distributed 
within the blazar sample. 

The lower peak frequency of FR-Is in vY^, could 
be understood as a result of different beaming fac- 
tors for the emission dominating in blazars and 
FR-I cores. This s eems to be in ag r eement with 
the model of e.g. IChiaberge et al.l (|2000bl ) who 



assume that the jet has a velocity structure with 
a fast moving central spine (which is strongly 
beamed and dominates in aligned sources) and 
a slower - but still relativistic - outer layer of 
more isotropic, less strongly beamed emission 
which dominates in the misaligned cases. Due 



'Using measurements for t he extended radio lum inosity as- 
sociated with the blazars l lAntonucci fc Ulvestad 1985) we 
introduce a luminosity threshold (and varying this thresh- 
old between log L1.4 GHz.ext of 31.5 and 32.5 erg/s/Hz) for 
identifying "FR-I blazars". Taking those blazars with ex- 
tened radio luminosities below the threshold as the possible 
aligned counterparts for our misaligned FR— I sources we see 
no trend of blazar FR type in Fig. llOl 
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Wills et al.„2004l 



to a shift in peak frequency of the emission, differ- 
ent amounts of relativistic beaming for aligned (i.e 
BL Lac) and non-aligned (i.e. FR-I) soures would 
also affect the observed broad-band spectral in- 
dices. This can serve as one possible explanation 
for the difference in a^o (or SE P curvature) be- 
tween FR-Is and BL Lacs (e.s. Chiaberge et al 




Trussoni et al. 2003: Balmaverde et al 



4.7. MIR and emission line luminosity 

Under the assumption that the narrow-line 
region of AGN is extended and not affected by 
nuclear obscuration it has been argued that the 
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Fig. 11. — MIR luminosity at 24 /iin plotted versus 
the [O III] emission- line luminosity. We sho'w the 
objects presented in iDicken et al. I (l2009l) as filled 
symbols and the FR-Is from this paper as open 
symbols. The FR-I points also indicate the domi- 
nating factor of the MIR emission: AGN - AGN- 
heated dust; synch ~ synchrotron emission; SF - 
star formation; host - stellar emission (plus pos- 
sibly star formation). Sources labeled "unkno'wn" 
(3C129 and 3C424) have MIR continuum emis- 
sion in excess of host galaxy contributions but due 
to missing CCC measurements we cannot securely 
distinguish bet'ween AGN or synchrotron emission. 
The dashed line is the bisector of a linear least 
square fit to the filled symbols (see text for de- 
tails). 

luminosity of the [O III] A5007 A emission line 
can serve as an isotropic tracer of the intrinsic 
radia t ive AGN luminositv (iRawlings fc Saunders! 

It has also been 
AGN the radiative 



1991: Tadhunter et al. 199: 



shown that for lumninous 
power, as traced by [O III] , is correlated with the 
MIR luminosity at 24 ^m |Tadhunter et al.l 12007 : 



In Fig.[TT]we plot the luminosity at 24 /im ver- 
sus the luminosity in the [O III] emission-l i ne for 
the the objects presented in iDicken et all (|2009l 
filled symbols) and for our FR-I sources (open 
symbols). We also show as a dashed line the 
bisector of a linear least square fit to the filled 
symbolf0. While the AGN FR-Is on average fall 
close to the bisector the visual inspection of Fig.fTTl 
also reveals that the bulk of the remaining FR-I 
sources tend to have lower 24 ^m luminosities. In 
fact, a Kolmogorov-Smirnoff test of the displace- 
ment in L24/im from the fitted line gives only a 
~ 1 % chance that two samples drawn from the 
same parent population would show a difference 
as strong as we see for the D09 sample and the 
non-AGN FR-I^ The FR-I AGN sources on 
the other hand are statistically consistent with 
being drawn from the same parent population as 
the D09 sample. We note that this comparison 
assumes that the 24 fim flux is dominated by an 
AGN component. While this is a reasonable as- 
sumption for the majority of the sources in the 
D09 sample (which mainly consists of NLRGs and 
quasars) we have shown here that it is certainly 
not the case for the non-AGN FR-Is. For the lat- 
ter the bulk of the MIR emission can be accounted 
for by other processes (Tab. [2]). 

The differences between FR-Is with and with- 
out warm dust can also be established by com- 
paring their respective MIR-to-[0 III] luminos- 
ity ratios: We find log [uL^ {24 ^im)/ L{[0 III])] = 
2.74 ± 0.06 for the four AGN FR-Is but only a ra- 
tio of 2.06 ± 0.20 for FR-Is without clear signs for 



di Sereeo Aliehieri et al] l ll997f) . ISimpsonI l ll998t) . See also 
Haas et al .1 l|2005l ) for an infrared/optical comparison. 



Dicken et al.1l2009H^ . 



^"The bisector was determined using o nly the filled s ymbol s 
and following the same strategy as ]Dicken et al,| j2009|) : 
Sources with upper limits on their [O III] emission were 
excluded as well as sources with clear starburst signatures 
in their optical spectra. For the rest of this section we refer 
to this sample of 33 objects as the "D09 sample". 

^^In order to match the selection criteria of those two samples 
we excluded FR— I sources with [O III] upper limits from 
the K— S statistics as well as sources which are dominated 
by star formation (plus BL Lac). 



^For distant FR-II radio galaxies and quasars an interest- 
ing thread in the literature offers good evidence for the 
anisotropy of [O III], which we do not pursue here since we 
concentrate o n FR-I objects. But thes e interesting refer - 
ences are e.g. I Jackson fc Brownd l|l990f ).lH es et al ] ||1993D . 
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Fig. 12. — Diagnostic op tical emission-line diagrams. On the left the dotted li ne is the divid i ng lin e used by 



KaufFmann et al. ( 20031 ) and the dashed line corresponds to the dividing line of Kewlev et alJ ( 200l[ ). Sources 



with [N II]/Ha>0.6 are often defined as Seyfcrts for [O III]/H/3>3 a nd as LINERs f or [O III]/H/3 < 3. For 
the diagrams in the middle and on the right the dividing lines from Kewlev et alJ ( 20061 ) are used. The 
symbols are the same as in Fig.[TTJ For 3C270 the data point with the higher [O III]/H/3 ratio corresponds 
to the small aperture HS T da,ta, while the lower flux ratio represents the ground-based results (see text). 



warm dust emissiorrn. The results of this exercise 
are just as we would expect from the spectral anal- 
ysis and provide a strong consistency check for our 
method. 

4.8. Diagnostic diagrams 

Optical emission-line diagnostics are com- 
monly used to identify the do minant activity 
type 



1981 



2001 



i n li ne-emitting sources (e.g. Baldwin et al 
IVeilleux fc Ostc rbrock 1987; Kewlev et al 



Kauffmann et al.ll2003l ). Although we here 
largely focus on the MIR continuum properties of 
our FR-I sources, a brief study of their optical 
emission-lines ratios with respect to their MIR 
spectral "type" seems promising. We take data 
for the emission-line ratios e ither from the refer- 
ences given in Ta b. [H or frompHansen et aD ( 19951 



3C218) and from lLewis et al.r(j2003l . IC4296). 

In Fig. [12] we show three emission-line ratio di- 
agrams for the FR-I objects (minus the broad- 
line sources 3C120 and E1821-I-643 and excluding 
BL Lac). Virtually all FR-Is, regardless of their 



^■^We here exclude the luminous broad-line AGN 3C120 and 
E1821+643, the synchrotron dominated sources BL Lac 
and 3C274 as well as the four sources which are dominated 
by emission from star formation (Tab. [2}. Note also that 
the slope of the fitted line shown in Fig. llll is not equal one 
meaning that the 24fiiai/[0 III] ratio changes with lumi- 
nosity. 



MIR properties, fall into the LINER or Seyfert 
region. Even the star-forming sources occupy 
the LINER region rather than showing Hll-type 
emission-line ratios (t his is discussed in the con- 
text of ULIRGs in e.g. lAntonucci|[2002al lbl). Only 
one source (3C386) occupies the transition region 
(Fig.ini left) or the HII region (Fig.[12l middle). 

For the FR-Is with warm dust emission, 3C15 
and NGC6251 fall into the Seyfert region for all 
three diagrams, which historically requires the 
presence of an active nucleus in these objects. 
3C84 and 3C270 on the other hand can be classi- 
fied as Seyferts or as LINERs depending on which 
diagram is used (Fig. fT^ . This, however, also 
shows that sources which are LINERs optically 
can have significant emission from warm dust in 
the MIR. Interestingly, both of these LINERs with 
warm dust have broad emission-lines detected: 



3C84 shows a broad Ha comp onent ()Ho et al 



Il997bl ) while iBarth et al.l (|l999f) report a (ten- 
tative) broad line in polarized ligh t for 3C270 



This i s consistent with the findings of lSturm et al 



( 20061 ) who show that type-1 LINERs on average 
have an additional warm dust component com- 
pared to type-2 LINERs. Therefore we can spec- 
ulate that LINERs with warm dust emission will 
also show broad emission lines. 

Because nuclear emission-line ratios can be di- 
luted by more extended emission from e.g. star 
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formation, especially for the weak and elusive 
sources considered here, we demonstrate the im- 
portance of aperture effects in the case of 3C270 
which appears twice on each diagram: For every 
line-ratio combination shown fFig.fT^ 3C270 is 
always classified as a LINER when using the the 
emission-line fluxes measured in a 2"x4" aper- 



ture (|Ho et al.lll997ar ). However, considering only 
the line ratios for an HST spectrurn of th e very 
nuclear 0.09" x 0.09" (jPerrarese et al.lll996l ) places 
3C270 in the Seyfert region for two out of the three 
diagrams. 

Fig. [12] also reveals that some of the FR-Is 
which are dominated by synchrotron emission in 
the MIR (and which lack clear signs for substantial 
warm dust) exhibit Seyfert-like emission-line ra- 
tios and many fall very close to the LINER/Seyfert 
transition line. 

5. Summary and conclusions 

The MIR spectra of 25 FR-I radio galaxies are 
presented. Many of these sources show contribu- 
tions from stellar emission at the shortest MIR 
wavelengths observed. The spectra generally turn 
over from a blue to a red continuum slope (in F^) 
around ^ 10 /im. Signs for star formation as traced 
by PAHs is detected in several objects and the rel- 
ative contribution to the total spectra range from 
not detected (e.g. 3C66B, 3C317) to minor (e.g. 
3C129, 3C270) to dominant (e.g. 3C31, 3C293). 

Focusing on sources with detected optical com- 
pact cores, the nuclear SEDs for 15 objects are fit- 
ted with a parabolic function in order to estimate 
the synchrotron contribution to the MIR emission. 
For 7/15 objects (47%) the nuclear MIR emis- 
sion turned out to be dominated by non-thermal 
processes: Apart from the stellar contributions 
at short wavelengths, the MIR spectra of these 
sources are consistent with the optical and ra- 
dio core measurements. Mid-infrared continuum 
emission in excess over the synchrotron cores is de- 
tected in the 8/15 remaining sources. For four out 
of those eight objects the spectra are dominated by 
processes related to the host galaxy (star forma- 
tion, stellar emission) . The remaining four sources 
(27%) - 3C15, 3C84, 3C270, and NGC6251 - 
show clear signatures of thermal MIR emission 
in excess over the non-thermal core contributions. 
This excess cannot be explained by either stellar 



emission or star formation, which argues for the 
presence of thermal dust emission powered by an 
AGN. At least one of these objects has been (ten- 
tatively) confirmed to have a hidden type-1 AGN 
with spectropolarimetry. 

Thus, for most of our sources in which the nu- 
clear MIR component can be identified (and is not 
masked by host galaxy emission) this emission is 
likely to be of non-thermal origin. But some FR- 
I sources do exist which show evidence for a hid- 
den AGN, consistent with findings at other wave- 
lengths. 

Parabolic functions were fitted to the radio and 
optical core data in order to estimate the nu- 
clear synchrotron component. These fits were per- 
formed in Fjy as well as in i^Fi, . Comparing the fit 
parameters for the FR-Is with those for a set of 
blazars and related objects, it appears that in F^, 
both types of objects show a similar relation be- 
tween the peak frequency and the curvature of the 
SED. However, while FR-Is cover about the same 
peak frequency range as the blazars used here 
for comparison, they show systematically stronger 
curvature in their SEDs. This reflects in the fit 
parameters for j/F^ where the FR-Is show on av- 
erage smaller peak frequencies and again stronger 
SED curvature. Different amounts of relativistic 
beaming in BL Lacs and FR~Is could serve as a 
possible explanation for these results. 

The comparison with a correlation between 
MIR and optical emission-line luminosity found 
for luminous AGN reveals that the AGN dom- 
inated FR-Is are consistent with that relation 
while the non-AGN FR-Is deviate substantially. 
In addition we show that FR-Is with and with- 
out warm dust differ significantly in their aver- 
age 24/zm/[0 III] luminosity ratio. The conclu- 
sion about the presence of an AGN in some of 
the FR-Is is further strengthened by their optical 
emission-line ratios. These tests provide a strong 
consistency check for our method. 

One important conclusion of this paper, as well 
as of other results in the literature, is that while 
most FR-Is lack powerful type-1 AGN, it is not 
tenable to generalize on associations between FR-I 
galaxies and "non-thermal only" AGN. A signifi- 
cant fraction of FR-Is do have warm dust emission 
which could by attributed to hidden type-1 nu- 
clei, a conclusion also a nticipated on other grounds 
(e.g. lAntonuccill2002a[ ). 
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A. Appendix: Examples for the subtraction of the star-forming template 

Here we present additional plots on the subtraction of the star-formation component. The object 3C270 
is already shown in Figs-O and [S] while 3C83.1 as well as 3C272.1 can be found in Fig. [SI This leaves 
3C129, 3C386, and IC4296 as objects in which a significant correction due to emission from star formation 
was applied. We do not show objects which are clearly dominated by star formation in their MIR spectra 
(see Tab. [5]). We also do not show objects where the contributions from star formation - as traced by the 
PAH features - is very minor and thus only a very small correction was made (3C84, 3C120, NGC6251, 
£1821+643; Tabs.[12]). 
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Fig. 13. — Subtraction of the star-forming template for 3C129 (left) and 3C386 {right). The top panel 
shows the observed spectrum with the scaled star-forming template. The bottom panel shows the corrected 
spectrum with a scaled spectrum of an early-type galaxy (NGC 1549). 
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Fig. 14. — Subtraction of the star-forming template for IC4296. The panels are organized the same way as 
in Fig.[T31 Here it becomes obvious what is mentioned in the note to Table [5] Due to the strong relative 
contribution from the synchrotron core (long-dashed line; Fig. [8]) the scaling of the templates to features in 
the spectrum can result in overcorrections. Also, the steep slope at A > 30 fim hints at substantial emission 
from cold dust (which is supported by FIR measurements; Fig. [7]) might influence the measurement at 30 /xm. 
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